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Abstract. The ^^C^^'N/^^C^^N and ^^C^^N/^^C^^N isotopic ratios have been determined in comets C/1995 Ol (Hale-Bopp), 
C/2001 Q4 (NEAT) and C/2003 K4 (LINEAR) at heliocentric distances of, respectively 2.7, 3.7 and 2.6 AU. These ratios have 
also been measured at r ~ 1 AU. No significant differences were found between all determinations, nor with the value obtained 
for other comets. If confirmed, the discrepancy between the nitrogen isotopic ratios from optical and millimeter measurements 
on CN and HCN would rule out HCN as a major parent of the cometary CN radicals. 
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1. Introduction 

■ Determination of the abundance ratios of the stable isotopes 
' of the light elements in different objects of the solar system 
provides important clues in the study of its origin and early 
history. Comets carry most valuable information regarding the 
material in the primitive solar nebula. 

We identified for the first time a number of emission fea- 
tures belonging to the ^^C^'^N B-X (0,0) band in spectra ob- 
tained with the SOFIN spectrograph at the Nordic Optical 
Telescope (Canary Islands). This allowed us to make the first 
optical measurement of the nitrogen isotopic ratio ^'^N/^'^N 
' in a comet (Arpigny et al. I2003t . This ratio was found to 
be lower by a factor of about two than the terrestrial value 
(272) and, less than half those obtained in Hale-Bopp from 
millimiter measurements of HCN, a possible main parent of 
CN (14n/15N= 323 ± 46 at r = 1.2 AU: Jewitt et al 119971 
330 ± 98, r = 0.92 AU: Ziurys et al. 1999 1. Spectra of fainter 
comets obtained with the ESO VLT and the 2.7m McDonald 
telescopes gave similar results, viz 140 ± 25 and 170 ± 50 for 
comets 122P/de Vico (r = 0.66 AU) and 153P/Ikeya-Zhang 
(r = 0.92 AU) respectively (Jehin et al. 2004i. On the other 
hand the optical determinations of the ^^C/^^C ratio consis- 
tently yield values around 90, only slightly lower than the HCN 
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Fig. 1. Spectra of comet C/2001 Q4 at 3.7 AU (top) and 0.98 
AU showing part of the R branch of the CN B-X (0,0) band. 
The corresponding quantum numbers are indicated at the top. 
The major differences are related to the heliocentric distance 
and show up in the excitation of lines of higher quantum num- 
ber when closer to the Sun. The intensity scale is relative only. 



millimiter measurements (respectively 111 ± 12 and 109 ± 22 
from the above-cited authors), and in agreement with the solar 
value (89). 

The differing results for nitrogen in CN and HCN would 
indicate that cometary CN radicals are produced from at least 
one other source with a much lower N isotopic ratio. 
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A monitoring of Hale-Bopp showed that somewhere be- 
yond about 3 AU, the production rates of HCN and CN agree 
with one another within certain limits (Rauer et al. pD03 ). 
Hence, sufficiently far away from the Sun, HCN could be the 
only source of CN such that ^^N/^'^N from CN should match 
the HCN value. In order to investigate this possibility, we mea- 
sured the N isotopic ratio in three comets at such heliocentric 
distances. By comparison, our previously published data con- 
cern comets much closer to the Sun (0.6-1.3 AU). 

2. Observations and data analysis 

Comets C/1995 Ol (Hale-Bopp), C/2001 Q4 (NEAT) and 
C/2003 K4 (LINEAR) have been observed at hehocentric dis- 
tances of, respectively, 2.7, 3.7 and 2.6 AU, and all three were 
also observed close to perihelion. Various telescopes were used 
in Chile, the USA, and the Canary Islands. Spectra of the 
CN violet band were obtained at resolving power of 60000 or 
higher (see TableQ. 

The techniques used to reduce the data and measure the 



isotopic ratios are described in Arpigny et al. (I2003> and Jehin 
et al. (2004, 1. The ^^c/i^c and ^^^N/^^N ratios are obtained si- 
multaneously by fitting the spectra with theoretical models. An 
empirical coiTection was adopted to take account of additional 
effects not covered by the pure fluorescence model (coUisional 
and possibly other unidentified effects). 

Indeed, the comparison of theoretical and observational 
^^C^*N data for several comets shows that a 3-parameter func- 
tion, with the quantum number, N, as the independant variable, 
always brings a satisfactory agreement. The intensity /jv of line 
N is corrected in the following way: 



F, 



N 



fNe^ 



1) 



(1) 



with a,b,c being estimated through an iterative process, for 
each individual spectrum. The same correction is applied to 
I3(^i4|y ^jj^ i2(-.i5]y mjjjgj- tjjg assumption that these spectra 
are similarly affected. These coiTections typically stay within 
15%. In cases of severe collision effects, in the central region 
of the coma of C/2001 Q4, they can reach up to 40%. One 
spectrum taken on May 6, 2004, right on the nucleus shows 
corrections of from —65% at R(0) up to +65% around R(15). 
The rms deviation of the fit between the observations and the 
coiTected model is generally between 5 and 15%, i.e., 3 to 4 
times smaller than for the uncorrected model. 

Only the R lines of the B-X (0,0) band (i.e., shortward of 
3875 A) are used since the P lines of the three isotopes are 
strongly blended. The spectra obtained at large heliocentric dis- 
tance have relatively low S/N ratios. This makes the use of a 
global fitting technique mandatory. As seen in Fig[T]the bulk of 
the CN emission is concentrated in fewer lines, namely those 
of low quantum number, when farther from the Sun. This in- 
creases the overall S/N ratio for the isotopic lines since the 
noise is dominated by the background and the CCD readout. 
Simultaneously, the low-excitation spectra are less contami- 
nated by the P lines of the (1,1) band of the dominant species. 

The estimated isotopic ratios are presented in Tableland 
plotted versus the heliocentric distance in Fig |2l The val- 
ues based on the Hale-Bopp SOFIN observations are shghtly 
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Fig. 2. Nitrogen and carbon isotopic ratios at various heliocen- 
tric distances 



different from (and supersede) those given by Arpigny et 
al. 



J2003t . This is mainly due to the use of better molecular 
parameters in the fluorescence model and of the ad-hoc empir- 
ical coiTection. The Hale-Bopp SOFIN values agree with the 
ratios from McDonald spectra obtained quasi simultaneously. 

3. Discussion 

Our measurements show that changes in the ^''N/^^N and 
^^C/^^C ratios with heliocentric distance are small as deter- 
mined from the CN radical. This means that the CN isotopic 
mixture remains constant within the explored range of helio- 
centric distances, and is the same for every comet, with average 
values of, respectively, 86 and 141 for ^^C/^^C and i^N/i^N. 
The carbon ratio is compatible with previous measurements, 
including millimiter measurements in HCN. On the other hand 
the nitrogen ratio is compatible with optical measurements of 
CN in other comets, but not with the HCN observations of 
Hale-Bopp. 

In order to get some idea of the constraints on the actual iso- 
topic variations, we assume two sources of CN, namely HCN 
and X. Let us call C = [i2ci4N/i2ci5N]totai the observed 
isotopic ratio of CN, and Ch = PC^N/i^ci^NIhcn the iso- 
topic ratio of the CN radicals coming from HCN. Similarly, we 
define Cx = [i^C^N/i^CiSNlx. 

We caU e = CNncN/CNtotai « i^CI^Nhcn/^'CI^N,, 
the fraction of CN coming from HCN. 



Hotal 



One has 

1/C«e(l/CH-1/Cx) + 1/Cx, 
so that 

(l/Ci-l/C2)«(ei-e2)(l/CH-l/Cx) 



(2) 



(3) 



where subscripts 1 and 2 refer to different observations — 
e.g., the same comet at different heliocentric or nucleocen- 
tric distances, or different comets. Hence, assuming that Cx 
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Table 1. Observational data, r is the heliocentric distance, is the heliocentric visual magnitude. 2dc is the McDonald 2D- 
coude spectrograph, SOFIN and UVES are the high-resolution spectrographs at, respectively, the NOT and ESO VLT telescopes. 
n is the number of spectra obtained, for a total exposure time of Exp. Extent is the size of the zone explored at the distance of the 
comet. R is the instrumental resolving power. Fwhm is the measured width of the CN lines. 





Date Comet 
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71lr 
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(km) 
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Fwhm 

(A) 


1 


Oct96 




2.73 


-20.7 


3.0 


2dc 


17 


30.6 


10000 


60000 


0.079 


2 


Apr97 


C/1995 Ol (Hale-Bopp) 


0.92 


3.0 


-1.3 
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2.8 


2000 


190000 


0.040 


3 


Apr97 




0.93 
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20000 


150000 


0.053:0.066 
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Sep03 


C/2001 Q4 (NEAT) 


3.73:3.66 


-18.8 


9.5 


UVES 
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23.4 


20000 
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May04 


0.98 


-5.4 ; -4.8 


5.6 


UVES 


11 


20.7 


50000 


80000 


0.052 


6 


May04 


C/2003 K4 (LINEAR) 


2.61 


-20.3 


9.2 


UVES 


1 


4.4 


7000 


70000 


0.056 






















7 


Nov04 


1.20 


14.7 


6.3 


UVES 


2 


3.0 


5000 


80000 


0.050 



and <^H do not vary, changes in the observed ratio l/( = 
i2c15N/12C^'*N are directly proportional to ei - £2, i.e., to 
changes in the relative importance of the column density of 
CNhcn or, with the opposite sign, of CNx. 

Let us consider an isotopic mixture with the observed ratio 
Ci = 140 and an isotopic ratio (n ~ 325 corresponding to the 
millimiter observations of HCN. Equation (|2j links (x and ei. 
Equation Q can be used to calculate the proportion £2 needed 
to produce a significant change in the observed ratio, say (2 = 
200, which would be easily detected. 

First, we assume that (i = 140 characterizes a mixture of 
CNhcn and CNx in equal proportion (£1 = 0.5). The isotopic 
ratio of X is then Cx ~ 89. The proportion of CNhcn must 
increase to £2 = 0.76 in order to yield an overall C2 = 200. 

As another (extreme) example, let us assume that the ob- 
served overall Ci ~ 140 characterizes CN coming exclusively 
from X (£1 = 0). The i^N/i^N ratio C2 ^ 200 is reached when 
53% of CN comes from HCN (£2 = 0.53). 

Our observations show that there is no such change in com- 
position up to 3.7 AU. 

Then, the hypothetical transition zone between a pure HCN 
source and a multiple source of CN must be farther out than 3.7 
AU. This does not fit well with Rauer et al.'s t2003j results. 

From Eq. (|3}, imposing the constancy of the isotopic ratios 
measured in several comets at various distances (Ci = C2), we 
find: 

(i) Cx = Ch = C (= 140), isotopically homogeneous sources, 
no constraints on the mixture proportions, or 

(ii) £1 = £2, the mixtures have identical proportions. 

We shall not discuss further hypothesis (ii) because a highly 
improbable conjunction of coincidences would be needed to 
achieve the constancy of the fractions e in every comet, under 
any conditions. 

The first hypothesis (i) implies that either the observed 
HCN or CN nitrogen isotopic ratio is faulty, or that the ob- 
served CN does not come from HCN. 



In the optical domain, the measurement of the CN isotopic 
ratios is straightforward as soon as the S/N ratio is high enough. 
The number of lines involved (a minimum of 7 in the present 
analysis) makes the possible influence of underlying unidenti- 
fied features rather improbable. The fluorescence models ap- 
pear to be reliable, so errors by a factor 2 seem to be excluded. 

The HCN isotopic ratio is difficult to measure, and so far 
has been obtained in only one comet, the exceptionally bright 
Hale-Bopp, but the two independent measurements agree with 
one another and there is no reason to discard them. If there 
was a problem, the carbon ratio should probably show it in 
similar proportions. Nevertheless, it would be important to see 
whether the HCN results can be confirmed on a broad sample 
of comets in order to pin-point the nature of the problem of the 
isotopic ratio of N. Are there actually large N isotopic differ- 
ences between the various constituents of a comet? Or does the 
cometary material have an overall excess of ^^N as compared 
to our cuiTent estimates of the solar nebula? Did the problem 
arise in this primitive nebula or is it linked to the formation and 
history of the cometary bodies? 

The discrepancy between the observed CN and HCN iso- 
topic ratios could be explained if most of the CN released by 
HCN is produced outside of the zones explored by the observa- 
tions, i.e., farther than about 5 10^ km from the nucleus. HCN 
would then photodissociate slowly, much farther than expected 
from the nucleus. Its contribution to CN along the line of sight 
toward more central regions of the coma would be negligible. A 
sufficiently low dissociation rate would make the correspond- 
ing CN invisible even at large heliocentric distances. However, 
the spatial distribution of HCN is rather steep and excludes 
such large scale lengths (see, e.g. Irvine et al. 1998 ). The major 
sources of CN — at small and large heliocentric distances — 
must then be some other molecule(s) with a dissociation scale 
length rather comparable to that of HCN and a lower ^"'N/^^N. 

Indeed, other sources of CN — simple or complex organic 
molecules — have been mentioned long ago in the literature, 
notably C2N2 (see, e.g.. Swings & Haser l 19561 . 
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Table 2. The carbon and nitrogen isotopic ratios, and the cor- 
responding estimated eiTors in parentheses. 



Comet 



C/1995 01 (Hale-Bopp) 



C/2001 Q4 (NEAT) 



C/2003 K4 (LINEAR) 



80 (20) 
90 (30) 
95(40) 

70 (30) 
90 (15) 

85 (20) 
90 (15) 



scale lengths and between production rates. This yields differ- 
ent results in different comets, a variability which would show 
up in the isotopic ratios. The isotopic measurements provide an 
independent and powerful fingerprinting method. If confirmed, 
they would show that HCN is not the unique — and probably 
not the major — parent of CN at small as well as large helio- 
centric distances. 

The origin of CN is still a puzzle, as is the apparent iso- 
topic differentiation within cometary constituents and its con- 
sequences on our understanding of the early solar system. 
Additional determinations of the N isotopic ratio in HCN and 
other N-bearing molecules are highly desirable in order to shed 
light on these issues. 



The difficulties met when investigating HCN and CH3CN 
as possible parents of CN were raised by Bockelee-Morvan & 
Crovisier ( 1985) on the basis of the expansion velocity in the 
coma. They suggested C2N2 and HC3N as possible alterna- 
tives. After re-examination of many near-UV observations of 
the CN radicals and high-dispersion spectra of the CN Violet 
(0-0) band in the coma of comet C/Austin 1989 XI, Festou 
et al. [T9981 favor a single parent which is neither HCN nor 
CH3CN, but a molecule having a lifetime of about 3.5 10'' s at 1 
AU. They propose the chemically stable C2N2 as a likely can- 
didate. Following the same lead, Bonev and Komitov J2000I I 
find that C2N2 is indeed to be prefeiTed over HCN as a sin- 
gle parent of CN. A recent analysis of the CN profile in comet 
21P/Giacobini-Zinnerby Lara et al. (500T) showed that HCN 
should be ruled out as the sole parent of CN. On the other hand, 
several analyses present HCN as the major or sole source of CN 
(e.g., Korsun & Jockers. l2002l Lara et al."2004'). 

Fray et al. (.2004) propose HCN polymers as a source of 
CN. A'Hearn et al. l2002| provide strong evidence that most CN 
are produced from grains in the coma rather than from nuclear 
ices. The link with CHON particles was suggested. However, 
the correlation between the dust and CN distribution in the co- 
mae is not perfect. This was shown by A'Hearn et al. (19^6) in 
the case of comet P/Halley. In a study of comet Hale-Bopp, 
Woodney et al. ( I2004> find that there is a better correlation 
between HCN and CN than between HCN and the optically 
dominant dust. A dust source could then be small or/and low- 
albedo undetectable particles instead of the observable dust, 
like black HCN polymers (e.g. Rettig, et al. I1992> . A very 
dusty comet like Hale-Bopp shows the same isotopic character- 
istics as quasi dust-free comets (e.g., 122P/de Vico, see Jehin et 
al. l2004l ) and the latter do not show abnormally low CN produc- 
tion rates. Greenberg and Li ( 1998 ) conclude that the source of 
CN and some other simple molecules is the organic component 
in comet dust. 

The IDPs with elevated ^^N point toward the presence 
of organic molecules, e.g., aromatic hydrocarbons (Keller et 
al. E004\. These may be related to the organic component of 
the cometary grains and explain the high ^^N abundance (see 
also Arpigny et al. l2003t . 

The link between CN and HCN in cometary comae has 
generally been investigated by looking for similarities between 
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